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We calculate the important quantity of superfiuid hydrodynamics, the relativistic entrainment 
matrix for a nucleon-hyperon mixture at arbitrary temperature. In the nonrelativistic limit this 
matrix is also termed the Andreev-Bashkin or mass-density matrix. Our results can be useful for 
modeling the pulsations of massive neutron stars with superfiuid nucleon-hyperon cores and for 
studies of the kinetic properties of superfiuid baryon matter. 
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^ I. INTRODUCTION 

Hi ! 

It is generally accepted that observations of pulsating neutron stars may potentially provide unique information on 
properties of superdense matter. For a correct interpretation of the observations one has to develop a hydrodynamic 
p_ _,, theory describing global pulsations. 

py^ In this paper we mainly focus on the hydrodynamics of superfiuid nucleon-hyperon matter of massive neutron-star 

| cores. More specifically, we study the important set of quantities of such hydrodynamics, the relativistic entrainment 
matrix. We assume the following matter composition: electrons (e), muons (/i), neutrons (n), protons (p), A-, and 
£~-hyperons (A and E, respectively). According to majority of calculations, A and E~ are the hyperons, which 
appear first in the stellar matter with the increasing density (see, however, Ref. [l[ and the discussion in Sec. I of Ref. 
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■£3 , Let us explain the physical meaning of the relativistic entrainment matrix. It is well known, that any baryon species 
i (i = n, p, A, or S) becomes superfiuid when temperature T falls below some critical value T C i. The microscopic 
calculations pit fil. IH la. I?! IH Itj, [To| predict that protons and hyperons pair in the spin singlet ^So) state, while neutrons 
in the core pair in the spin triplet ( 3 P2) state. In the most general case, when all baryon species are superfiuid, we 
have four condensates and a normal (nonsuperfluid) liquid component. The normal component includes electrons, 
muons, and baryon Bogoliubov excitations. We assume that, because of collisions, all these 'normal' particles have 
the same hydrodynamic velocity. Below in this paper we always work in the frame, comoving with the normal liquid 
, component. 

If there are no superfiuid currents in the system, then Cooper pairs are formed by particles with strictly opposite 

^ I \ momenta (e.g., p and — p). The presence of superfiuid currents corresponds to a situation when pairing occurs between 

particles with momenta (p + and (— p + QJ. The total momentum of a Cooper pair is then 2Q i . In this case the 
0> ' particle current density j i of a species i can be written as 0, [H| 

o ' 



> 
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c 2 J2^Q k - (i) 



Here the summation is performed over all baryon species, fc = n, p, A, and S; c is the speed of light. Finally, Yg. is the 



4x4 relativistic entrainment matrix. It is symmetric, Yg. = Y^, and generally depends on the four baryon number 
densities rii and temperature T. It follows from Eq. ([T]) that the superfiuid motion of one species contributes to 
particle current density of another species (and vice versa). This 'entrainment' effect was first suggested by Andreev 
and Bashkin [jjj in the context of superfiuid solutions of 3 He in 4 He. 

In the nonrelativistic limit the matrix Y^ (more precisely, its analog, the matrix pn~) was anal yze d for a neutron- 
proton mixture at zero temperature in Refs. [l3l . [l4l . [l5| and at arbitrary temperature in Ref. [16||. In our recent 
paper Q (hereafter GKH09) we calculate the relativistic entrainment matrix for the nucleon-hyperon matter in 
the approximation of zero temperature. The aim of the present study is to extend the analysis of GKH09 to finite 
temperatures. For that, the strongly interacting nucleon-hyperon mixture is considered in the frame of relativistic 
Landau theory of Fermi liquids, generalized to allow for possible superfluidity of baryons. 

In Eq. |T]) it is assumed that the matrix elements are scalars. Strictly speaking, this is the case only if the 
matter is isotropic in the absence of superfiuid currents, i.e. if all baryon species pair in the spin singlet ( x So) 

state. Meanwhile, neutrons undergo triplet-state pairing. This leads to anisotropic neutron energy gap Ap"' and to 
appearance of the preferred direction along the neutron quantization axis. As a consequence, the matrix elements Yik 
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are tensors with respect to spatial rotations (rather than scalars, as in the isotropic case). This makes the problem of 
calculation of the matrix Y ik much more complex and model-dependent [la ] . 

To avoid this difficulty we assume, following Refs. [IB, 53]' that the unperturbed matter can be treated as a collection 
of microscopic domains with chaotically directed neutron quantization axis. After averaging over the large number of 
such domains, the elements will become scalars. One can obtain then qualitatively correct results for the averaged 
matrix elements Yik, assuming pairing of neutrons, and introducing an effective isotropic energy gap A^g (T) in 
the neutron dispersion relation, 

AW(T)=min{A(")(|p|=p Fn )}. (2) 

Here A^(T) is defined as the minimum of the angle-dependent gap Ap on the neutron Fermi surface, see Refs. 
[l6l \vf\ for more details. Below in this paper we follow this strategy. 

The paper is organized as follows. In Sec. II we formulate a relativistic Hamiltonian describing superfluid nucleon- 
hyperon mixture. In Sec. Ill we diagonalize this Hamiltonian, determine baryon dispersion relations in the presence 
of superfluid currents, and calculate the relativistic entrainment matrix Y^. Section IV presents summary. 

Throughout the paper, unless otherwise stated, indices i and k refer to baryons, i, k — n, p, A, and S. We use 
the system of units in which the Planck constant ft, the speed of light c and the Boltzmann constant /cb equal unity, 
h = c = Are = 1. 



II. THE CHOICE OF A HAMILTONIAN 



Let us obtain a model Hamiltonian for a homogeneous system of relativistic, degenerate, strongly interacting baryons 
with superfluid currents. A natural way to obtain such a Hamiltonian is to use the framework of Landau Fermi-liquid 
theory, generalized to account for superfluidity by Larkin and Migdal [IH and Leggett • Following Ref. [n|, the 
Hamiltonian H describing a superfluid nucleon-hyperon mixture can be generally presented in the form 

H — ^iNi = HlF + Hp a i r i n g . (3) 

i 

Here, Nj and fii are, respectively, the number density operator and the chemical potential for baryon species i = n, 
p, A, E; Hlf is the Fermi-liquid Hamiltonian for the mixture; H pa i r i ng is the pairing Hamiltonian. 
If our system were nonrelativistic, Hlf would be given by (see, e.g., [HI, fl3|) 

Hlf = E £ } (P) (4 4)t 4 4) -^ } ) 

psi 

+ \ E f ik M) (W - 0? )f # - ( 4 ) 

pp' ss'ik 

Here p and p' are the particle momenta; s and s' are the spin indices; a p ^ = dpi = or aty is the annihilation 
operator of a Landau quasiparticle (not the Bogoliubov excitation!) of a species i in a state (ps). We restrict 
ourselves to a spin-unpolarized nucleon matter. This allows us to simplify the notations by suppressing the spin 
indices, whenever possible. Furthermore, 9 p ^ = (p Fi — \p\), where 9{x) is the step function; (p) = v Fi (|f>| — Ppi), 
where v Fi = p-pil m * and are i respectively, the Fermi- velocity and Fermi-momentum with m* being the effective 
mass. Finally, f lk (p,p') is the spin-averaged Landau quasiparticle interaction (we disregard the spin-dependence of 
this interaction since it does not affect our results). In the vicinity of the Fermi surface the arguments of the function 
f lk (p,p') can be approximately put equal to p « p M and p 1 « p Fk , while the function itself can be expanded into 
Legendre polynomials Pi(cos9), 

i 

where 9 is the angle between p and p' and f\ k are the (symmetric) Landau parameters, ff* = f kl . 

The Landau theory was extended to relativistic Fermi- liquids by Baym and Chin [2(|. These authors showed 
that the formal structure of the theory remains practically the same as in the nonrelativistic case. In particular, 
the expression for the variation of the energy density has the same form. Thus, we assume that the Hamiltonian 
describing the relativistic mixture is still given by Eq. We emphasize that this is our assumption; it would be 
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interesting to derive this Hamiltonian from the microscopic theory in analogy to what was done by Leggett [19( for 
the case of nonrelativistic Fermi-liquid. 

The relativistic effects affect only those properties of the theory which are related to transformation of various 
quantities from one frame to another. For example, as it is demonstrated in Ref. 20], the transformation law for 
the Landau quasiparticle interaction f lk (p,p') is more complex than that in the nonrelativistic theory. Also, in the 
relativistic theory the expression for the effective mass m* as a function of Landau parameters f[ k should be modified 



am, 



Mi 

m* 



= 1 - 



k 



(0) 



To reproduce the nonrelativistic result, one has to replace everywhere in this equation ^ with the mass of a free 
particle mj. In Eq. ([6]) m is the number density of particle species i, 



3tt 2 



(7) 



and the symmetric matrix Gik is defined by 



G, 



2 fik 
g n 4 ^Fi^Fk J 1 



(8) 



Let us obtain now an expression for the relativistic Hamiltonian H pa i r i ng , describing pairing between the Landau 
quasiparticles. The pairing in relativistic systems was studied in the literature mainly with application to quark 
matter (see, e.g., the early review [2l[ and a recent review [22| and references therein). However, there are also 
several papers that explore the importance of relativistic effects in terrestrial superconducting materials (see, e.g., 

Refs. Ha [Mm she mm, 

We start with the analysis of a pairing Hamiltonian for a mixture of relativistic noninteracting baryons. In other 
words, we ne glec t for a moment the Landau quasiparticle interaction f tk (p,p') in Eq. (0| for Hlf- As it is demonstrated 
in Refs. (23. 124, 12H [29j | . in the mean-field approximation the relativistic analog of the well-known non- relativistic 
Hamiltonian, responsible for 1 Sq pairing of particles, has the form 



H 



pairing 



\ E/ d3rdV AW (r,r') 



* W (r) 7 



5 v^V)-I <|* (i) (r) 7 5 4V)|)]+H.c. 



(9) 



Here A^(r,r') = A'''(r',r) is the superfluid order parameter; \&W(r) is the relativistic fermion field operator and 



*«(r) 
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¥°(r) 



is the charge conjugate field (see, e.g., the textbook [3l| for details). Furthermore, W (r) = 
\Er(*)t( r )^,o. ^5 _ ^0^,1^,2^,3^ w ]j ere y are the Dirac matrices (I — 0, 1, 2, and 3). The second term in square brackets 
in Eq. @ is presented to avoid double counting in the expression for the system energy. Terms of this kind is a typical 
feature of the mean-field approach. Similar terms appear in the mean-field formulation of the nonrelativistic pairing 
Hamiltonian (see, e.g., a review [321]). 

The Hamiltonian ([9]) can be simplified if we take into account that for a homogeneous system with superfluid 
currents the order parameter A^(r,r') can be written in the form (see, e.g., [2l| ) 



A«(r 



i) e ip(r-r') e iQi(r+r') 



(10) 



Here Ap is the Fourier component of the order parameter. It is even in p, 

A« = A (i) 



(11) 



since A^ (r, r') is a symmetric function of r and r' . Furthermore, 2Q i is the momentum of a Cooper pair in condensate. 
Below in this paper we always assume that Qi/p^k ^ 1 anc l restrict ourselves to a linear approximation in Qi/p Fk - 
Substituting Eq. (fTO)) into Eq. ([9]), one gets, after some algebra [with an accuracy up to quadratic terms in Qi/p Fk ], 



H, 



pi 



Ji)1 



-p+Q,l 



(I a, 



(Ot 



,(0t 



'P+Ott -p+Q,l 



„(<)t 



(0t J<)t 

'p+QiT -P+Q4 



H.c. 



(12) 
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The analogous derivation of a Hamiltonian for a homogeneous system is explained in detail in Ref. [28| for the case of 
Q i — 0. In Eq. (|12p Cp ^ is the creation operator for antibaryons of a species i. Generally, because H pa j r i ng does not 
contain 'interference' terms (e.g, of the form a p +Q j C p+Q .j_ ) > one ma y diagonalize the Hamiltonian Hlf + H pair i ng by 
simply performing the standard Bogoliubov transformation separately for baryons and antibaryons (this procedure is 
very clearly described in Ref. [28|; in this case the species index i in Eq. for Hlf should also run over antibaryons). 
However, since the conditions in the neutron-star cores are such that the population of antibaryons is negligible, one 
can neglect them and rewrite H pa i r j ng in its final form 



Hpairing — ^p 



i) 

pi 



„Wt „Wt 1 /|_W1 -001 II ] i n m \ 



We see that the Hamiltonian Hp airing for a homogeneous system of noninteracting relativistic baryons is given by the 
same equation as in the nonrelativistic limit. Moreover, it coincides [l6l with the pairing Hamiltonian describing 
strongly interacting nonrelativistic Fermi-liquid [the operators in Eq. (|13|) refer then to Landau quasiparticles rather 
than to real particles]. Thus, it is reasonable to assume, that H pairing for interacting relativistic baryon mixture is 
also given by Eq. (fT3"]) . All calculations below are made under such assumption. 

In the Appendix we give further arguments supporting the form of the expression (fl"3]) for H pa i r i ng . Namely, we apply 
the consideration given in this section to a specific case of relativistic mean-field model in which baryon interactions 
are mediated by various types of meson fields. We show that for this model the homogeneous Hamiltonian (|13p can 
be directly calculated from the general expression ^j. 

In view of Eq. pip , in the linear approximation the superfluid order parameter A p l \ entering Eq. (|13p . is the same 
function of p as in the system without superfluid currents. It can be chosen real, just as for the system with Q i = 
(e.g., [HI]). In this case Ap ^ is the energy gap in the dispersion relation for Bogoliubov quasiparticles (see below). 



III. THE RELATIVISTIC ENTRAINMENT MATRIX AT FINITE T 



A. General equations 

It follows from Eqs. (jlj and (fT3|) that the relativistic expressions for Hlf and H pa i r i ng are essentially the same as 
in the nonrelativistic case. Correspondingly, the further consideration is quite similar to that in Ref. [161 ]. Thus, here 

we only briefly sketch the main results, referring the reader to Ref. [lq for more details. 

(i) 

Introducing the Bogoliubov operators b P s 

n {l) - v (t) b il) +7; (l V i)f (U) 
p+Qil ~ U P °p+Q,]^ v p °-P+QiV { ' 

p+Qil ~ P °P+Q,l V P °-p+Q,T' ^ 10 J 

where u p ^ and v p ^ are even functions of p related by the normalization condition, 

u^ + v^ = l, (16) 

one obtains for the system energy density E 



'p+Q, °p+Q, 



psi 



+1 £ ^(p+Q^+ft)^-^,)^ 



(fe) _ Jk) 

+Q k V+Q» 



pp' ss' ik 



pi 



Here Np+Q. an( l -^p+Q arc ^ ne distribution functions for Landau quasiparticles and Bogoliubov excitations with 
momentum (p + Qj), respectively, 



AT« _ (I (')t (i) |\ _ (I (')t Ji) 



P+Qi - M "p+Qj "p+Qj I ~ ^P+QdS+Qd 



= ^ )2 + u« 2 4% i -,« 2 ^ p+0j , (18) 
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Notice that, since the parameters Up 1 and v p ^ are even functions of p, they do not depend on Q i in the linear 
approximation. 

The entropy of the system is given by the standard combinatorial expression, 



S=- 



E K 1 - ^U) ln (l - + ^~p+Qi ■ (20) 



pst 



Using Eqs. (|16 j) — (|20 jt and minimizing the thermodynamic potential F = E — \iirii — TS with respect to and 



CO 

Up ' , one gets 



■Fp+Q, ~ e m it > ( 21 ) 
1 + e Vk>»/ j 



1 / ffW _i_ rr(») 

Here <£ p +Q is the energy of a Bogoliubov excitation of a species i with momentum p + Q i ; Hp+Q is the energy of a 
Landau quasiparticle in normal (nonsupcrfluid) matter, 

H pl Qi = 4° (P + Q<) + E (P + Qi P' + Qk) «3-Q„ ^'IqJ ■ ( 24 ) 

p's'k 

Since Q j <§C p Ffc , it can be expanded in powers of Q^, 

H^ Qi =s^(p) + AH^\ (25) 

where e^(p) sw « K ([p| — p Fi ) is the quasiparticle energy in the absence of superfluid currents and AHp^ is a small 
current-dependent term. In the linear approximation AHp can be generally written as 

AiJ« = 7*fc pQ k , (26) 
fc 

where 7^. is a 4 x 4 matrix to be determined below (notice, that the definition of 7^ differs by a factor of rrik from that 

used in Ref. [IB]). Substituting Eqs. f|22|) and (|25[) into Eq. (|23[) one verifies that Up (and Vp^) is indeed independent 
of Q k in the linear approximation. 

The distribution function Np+Q. can also be expanded in powers of Q { . Using Eqs. (fT8| and ([2"l" |) -(f23 |) . one obtains 

^Sq^^ + ^jA^). (27) 

Here £p 1 ', fp\ and rip' denote, respectively, the quantities £ p +Q.> -^p+Q.) an d ^p+Q.' m tne absence of superfluid 
currents. They are given by the following well-known expressions: 



Ef = ^/ £ W2(p) + AW 2 , (28) 

fp° = ^, (29) 

4 4) = ^ )2 + (^ )2 -^ )2 ) (30) 
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Calculation of the matrix 7^ 



Let us calculate the matrix 7^. Substituting Eqs. (|25j) and (|27|) into Eq. (|24l) . one obtains, with the accuracy to 
the terms linear in Q { , 



2 p's'k \ 



9 ^) A ^' (31) 



The functions in curly brackets have sharp maximum near the Fermi surface of particle species k (at p' ~ p Ffc ) 
that the integrals in Eq. (f3"Tj) can be approximately calculated and presented in the form 

f (fe) 



E/*(P.JO^-0* = ~PQu. (33) 



£ /* (P,ff) r% A^ fe) = AffW , (32) 



dp' n, 

p's' ^ 

Here the matrix Gik is defined in Eq. ([8]) while the function <3>i is given by 

(34) 

It changes from <F; = at T — to <&i = 1 at T > T C i (we remind that T C i is the critical temperature for transition 
of particle species i to the superfluid state) . Assuming that the energy gap A^' does not depend on the momentum 
p, this function was calculated numerically and fitted in Ref. [34j ]. The fitting formula and further details concerning 
$i can also be found in Ref. (l6| . 

Using Eqs. ([26)) . (|32| . and p3|) and equating prefactors in front of the same Qfc in Eq. (|3ip. one obtains the following 
system of 4 x 4 linear equations for the matrix 7^ : 

&ik Gik v-* Gu * t . oc x 
7ife = — + — — m, * l7ll , (35) 

Notice that this matrix is not symmetric, 7^ 7^ 7^. In Eq. (|35[) for each fc we have four equations (with i = n, p, 
A, and E). These equations can be decoupled from the whole system and solved separately. Thus, there are in fact 
four independent systems of four linear equations. Their solutions can be immediately found in a number of limiting 

cases: 

(z) For very low temperatures, when all = 0, 

$ik Gik tne\ 

Hk = — + ■ (36) 

(ii) For nonsuperfluid matter, when all $j = 1, 

Hk = %• (37) 

m , 

(in) If the species in the mixture do not interact with each other, that is ff° — for i 7^ k, one has [see Eqs. © 
and ©] 

n i &ik /ool 

fit (rii + G u m* $0 

(iv) If $i = for all species except for the species j (i.e. $j 7^ 0), then 

Sik . G ifc GijiGjkrrf+njSjk)®, 
7tfc = — r H ? T7^ T^n — ■ (39) 

771* 77, 77j (72j + Gjj 777* $j) 
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(v) Finally, for any nonsuperfluid species j (for which = 1), the matrix elements 7^ related to this species, are 
given by the same Eq. (|37[) as in the case of a completely nonsuperfluid mixture, 

lij = % (40) 

TO* 

The general solution for the matrix ji k in case of four particle species is rather lengthy, but can be easily obtained 
from Eq. (|35[) and, for example, incorporated into a computer code. Here we present the solution for an important 
case of a mixture of two particle species (e.g., neutrons and protons) 

{n t + Gu to*) (n fc + G kk m* k $ fc ) - g| m* to* $ fc 
% = , (41) 

7. = S*"^ 1 -**) , (42) 



where 



5' 



5 = (n. t + G lt to* $0 (n fe + G kk m* k $ fc ) - G? fe to* m£ $ t $ fc . (43) 



In Eqs. (l4Tj)-([43l) indices i and k refer to different particle species, i ^ k. The presented solution agrees with that 
given in Ref. jl6l | for the nonrelativistic case (we remind the reader that the definition of 7^ in Ref. [l(| slightly 
differs from ours, see the comment to Eq. (|26p ). However, in contrast to Ref. [l6j . in Eqs. pTj) - P3")) one should use 
the relativistic expression for the effective masses given by Eq. ([6]). 



C. Calculation of the matrix Yik 



To find the relativistic entrainment matrix Y ik let us calculate the particle current density j t . For that, we make 
use of the fact that, as was emphasized by Le ggett 11911 . the expression for j i has the same form as in the case of usual 
nonsuperfluid Fermi- liquid (see also Refs. [2[ Il3l. Il4l. Il5l. [l6j ) . In other words, 

ps 

Substituting expressions (|25|) and (|2T[) into this equation and performing a simple integration, one obtains Eq. ([1]) 
with 

Y ik =n ilik {l-$ i ). (45) 

One may check that this matrix is symmetric, Y ik — Y ki . At zero temperature this equation can be simplified with 
the help of Eq. (f3"fi|) , so that one reproduces the result of GKH09, 

Y ik = ^6 ik +G ik . (46) 

TO 4 

Notice that, as follows from Eq. ©, in this limit the matrix Yi k satisfies the sum rule (see GKH09), 

^2^kY ik = m. (47) 

k 

In the nonrelativistic limit the matrix Y\ k is related to the nonrelativistic entrainment matrix pi k by the equation 
(see, e.g., [ll|) 

p ik = mim k Y ik . (48) 
Using this equation and Eq. (|45[) . one may reproduce the various limiting formulas for pi k , presented in the literature 

[SOS [Hi 

Our results are illustrated in the Fig. 1. The normalized symmetric matrix Yi k /Y is shown as a function of 
temperature T for the baryon number density rib = 3no = 0.48 fm -3 . To plot the figure we employed the third 
equation of state of Glendenning [35j | . The Landau parameters /f of nucleon-hyperon matter were calculated for this 
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log 10 T (K) log 10 T (K) 

FIG. 1: Normalized symmetric matrix Yit/Y as a function of T for the third equation of state of Glendenning [35|. The 
normalization constant Y = 3no//i n (3no) = 2.48 x 10 41 erg -1 cm -3 . Solid lines show the elements of the matrix Yik/Y; each 
curve is marked by the corresponding symbol ik [i, k — n, p, A, E). Vertical dotted lines indicate baryon critical temperatures. 

equation of state in GKH09. The normalization constant Y is taken to be Y = 3i%o/[i n (3no) = 2.48 x 10 41 erg -1 
cm -3 , where no = 0.16 fm -3 is the normal nuclear density and /x n (3no) = 1.94 x 10 -3 erg is the neutron chemical 
potential at n& = 3no- We choose the baryon critical temperatures T C i (i — n, p, A, S) equal to: T cn = 5 x 10 8 K, 
T cp = 2x 10 9 K, T cA = 3 x 10 9 K, and T cS = 6 x 10 9 K. 

On the left panel we plot the 'nucleon' matrix elements; the 'hyperon' matrix elements are shown on the right 
panel. In agreement with Eq. (|45p. the matrix elements, involving a nonsuperfluid species, vanish. For instance, if 
neutrons are normal (i.e., T > T cn ) 1 then $„ = 1 and Y n i — Yi n = 0. At T < 10 s K all the matrix elements approach 
their asymptotic zero-temperature values [see Eq. ([46]) ]. As follows from the figure, the nondiagonal matrix elements 
are essentially smaller than the diagonal ones and are all negative for the chosen equation of state. 

IV. SUMMARY 

The present paper is a continuation of GKH09, where the relativistic entrainment matrix Y^ of nucleon-hyperon 
mixture was calculated for the case of zero temperature. Here we extend the results of GKH09 to finite temperatures. 
For that, we employ the relativistic Landau Fermi- liquid theory [2(|, generalized to allow for superfluidity of baryons 

[SGI. 

We demonstrate that, as in the case of nonrelativistic neutron-proton mixture [l6[, the matrix Y^ is expressed 
through the Landau parameters f{ k and the function of temperature <&i. The Landau parameters f\ k for the relativistic 
nucleon-hyperon matter were derived in GKH09. The quantity <&i(T) is the universal function of T under assumption 

that the baryon e nerg y gap Ap is independent of momentum p. In the latter case the fitting formula for can be 
found in Refs. fillip. 

Our consideration of Y^ differs from the previous calculations of the entrainment matrix, available in the neutron- 
star literature [H, [HI, HE , by the following: 

(i) we consider the problem in a fully relativistic framework; 

(ii) we allow for the presence of two hyperon species (A- and S~-hyperons), in addition to neutrons and protons. 
It should be noted that our results can be easily extended to describe any number of baryon species (not necessarily 

four). The main problem then will be to determine the Landau parameters for these species. 

The calculated relativistic entrainment matrix is an essential ingredient in hydrodynamics of superfluid mixtures 
[36j . It can be important for studying the pulsations of massive neutron stars with superfluid nucleon-hyperon cores 
(see Ref. [37j for an example of such study). Also, since the matrix enters the expression (|2T|) for the equilibrium 
distribution function of Bogoliubov excitations, it can influence various kinetic properties of superfluid baryon matter, 
for example, the shear viscosity. The related problems will be analyzed in a separate publication. 
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Appendix 

The general approach of Sec. II can be illustrated with a specific example of a relativistic mean-field model in which 
interactions between baryons are mediated by meson fields. For definiteness, following GKH09, we consider a-w-p 
version of the mean-field model with self-interactions of scalar er-field. However, our consideration remains essentially 
unaffected if one allows for additional meson fields (e.g., <5-meson or hidden strangeness a* and </>-meson fields). 

In GKH09 wc reformulated the a-uj-p mean-field model in terms of the relativistic Landau Fermi-liquid theory. In 
particular, we calculated the Landau quasiparticle interaction function f lk (p,p')- Thus, we determined the Hamil- 
tonian Hlf for this model [see Eq. Q], Let us now turn to the Hamiltonian H pa i r i ng . It is still given by Eq. © 
with 

*W (r) = "i { a Ps « W fa s ) ^ + 4^ " W (P. s ) e " Pr ) . ( 49 ) 

ps * 1 

*W(r) = E 7^ ( a Ps 1 " W fa s ) c ~ 4pr + $ « W fa s ) eipT ) ■ ( 5 °) 

ps * 1 

The bispinors it^ (p, s) and i/W (p, s) correspond, respectively, to particles and antiparticles. For the mean-field model 
they can be written out explicitly (see, e.g., Ref. [381]). 



u w (»,«)= / ; ; , v w (p,s) = [ v v _ , - . (51) 

V^" M * (^W V + M t w -s 



In Eqs. (j49]) ~ ([5T|) £i — y P + M* 2 ; unit vector n is directed along P. a = (cr x ,a y ,a z ) is the vector composed of Pauli 
matrices; w s and w s are the spinors which are defined as in Ref. [39| (see §23 of this reference). Furthermore, 

P = P - guiv - 9 P ihiPz, (52) 
M* = mi- g al a, (53) 

where u>, p 3 , and a are the meson fields that are generated by baryon currents and densities; g U i, g P i, and g a i are the 
coupling constants; and 2^ is the isospin projection for baryon species i (for more details, see, e.g., GKH09). 

Substituting now field operators (pi9"j) and (j5U)) into the expression ^ for H pair i n g, one verifies, that the homogeneous 
pairing Hamiltonian for the mean-field a-uj-p model is given by the same expression as was obtained in Sec. II [see Eq. 
(TT2")) of this section or Eq. (|13p , if we neglect antibaryons] . This particular example presents an additional argument 
supporting the conjecture that the general expression for H pa i ring in the relativistic Landau Fermi-liquid theory, is 
also of the form (fT3"|) . 
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